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The main obstacle in describing inclusive decay spectra in QCD — which, in particular, limits 
the precision in extrapolating the measured B — > X s "/ rate to the full phase space as well as 
in extracting \V u b\ from inclusive measurements of charmless semileptonic decays — is their 
sensitivity to the non-perturbative momentum distribution of the heavy quark in the meson. 
We show that, despite this sensitivity, resummed perturbation theory has high predictive 
power. Conventional Sudakov-resummed perturbation theory describing the decay of an on- 
shell heavy quark yields a divergent expansion. Detailed understanding of this divergence 
in terms of infrared renormalons has paved the way for making quantitative predictions. In 
particular, the leading renormalon ambiguity cancels out between the Sudakov factor and 
the quark pole mass. This cancellation requires renormalon resummation but involves no 
non-perturbative information. Additional effects due to the Fermi motion of the quark in 
the meson can be systematically taken into account through power corrections, which are only 
important near the physical endpoint. This way the moments of the B — > X s ~/ spectrum with 
experimentally-accessible cuts — which had been so far just parametrized — were recently 
computed by perturbative means. At Moriond these predictions were confronted with new 
data from BaBar. 



1 Introduction 

Experimental measurements of inclusive decays, such as B — > X s ^f and B — ► X u lP, cannot 
be performed over the whole phase space. Therefore, in order to make precision tests of the 
Standard Model, we must have a quantitative understanding of the spectra. 

Most of the rate of these decays is associated with jet-like kinematics, where the invariant 
mass of the hadronic system is small. It is well known that this important region is difficult to 



describe theoretically [1—3]: the observed spectrum is directly sensitive to soft radiation off the 
heavy quark which is hard to disentangle from the "primordial" momentum distribution of the 
heavy quark in the meson. Technically, this sensitivity leads to the breakdown of the perturbative 
expansion in powers of a s (m b ) as well as that of the local Operator Product Expansion (OPE) 
in powers of A/m b . 

Let us consider the photon-energy spectrum in B — > X s j. The non-perturbative nature of 
this distribution is obvious from the fact that the physical endpoint is near E^ = Mb/2, where 
Mb is the meson mass, while at any order in perturbation theory — where the initial state is 
represented by an on-shell heavy quark — the spectrum vanishes for EL > m b /2, where m b is the 
quark pole mass. In reality the quark is off shell and its Fermi motion translates into smearing 
of the spectrum up to E^ ~ Mb/ 2. Given this physical picture the perturbative approach 
appears inappropriate. Consequently, the predominant strategy has been to parametrize the 
spectrum based on available data, rather than to compute it. On the other hand, since total 
inclusive rates, as well as their first few moments, are theoretically described by the OPE [4,5] - 
and in particular, if 0(A 2 /m|) terms are neglected, by perturbation theory — one expects just 
small, computable deviations when moderate kinematic cuts are applied [6]. Progressing to more 
stringent cuts, resummed perturbation theory is needed. Eventually it must be complemented by 
non-perturbative corrections. We will show that the predictive power of resummed perturbation 
theory is significantly higher than it superficially appears to be. 



2 The divergent nature of the perturbative expansion of spectral moments 

A key ingredient in the perturbative approach to inclusive spectra is the use of moment space. 
Considering the B — ► X s -y photon-energy spectrum in terms of x = 2E^/m b , the perturbative 
coefficients have singular structure at the perturbative endpoint, x = 1. Owing to cancellation 
of logarithmic infrared singularities between real-emission and virtual corrections, the moments 
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are infrared- and collinear-safe. has a well defined perturbative expansion to all orders. This 
expansion, however, does not converge well. It has: (a) Sudakov logarithms, lnA^, that are 
the finite reminder in the cancellation of logarithmic infrared singularities. These logarithms 
appear at any order in the expansion and dominate the coefficients at large N (b) infrared 
renormalons that reflect power-like sensitivity to soft momenta and induce factorial growth of 
the coefficients at large orders that renders the perturbative expansion non-summable. 

To make quantitative predictions both these sources of large corrections need to be re- 
summed. The resummation of Sudakov logarithms is possible thanks to the factorization prop- 
erty of matrix elements in the soft and the collinear limits. In decay spectra there are two 
distinct regions of phase space giving rise to infrared logarithms [7]: soft radiation with mo- 
menta of order m/N from the nearly on-shell heavy quark and fragmentation of the final-state 
jet of mass squared m 2 /N. In moment space these contribution factorize: 

= H(m b )J(m 2 h /N- M )S PT (m 6 /iV; fj) = H(m b ) x Sud(A^, m b ), (2) 

up to corrections 0(1/N). Here fj, is a logarithmic factorization scale. The Sudakov factor, 
which sums up the logs to all orders in a s (m b ), takes the form of an exponential: 
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"This smearing is related to the structure of the initial state, and it is independent of the appearance of 
resonances in the final state. 



where A = (/3 af s (ml) /it) In IV. Note that Sud(N,m b ) depends on the quark mass only through 
the argument of the coupling and that it does not depend on any factorization scale. The first 
three towers of logarithms in the exponent, i.e. c ntTl+ i — leading logs (LL), c n ^ n — next-to- 
leading logs (NLL) and c nin _i — next-to-next-to-leading logs (NNLL) are known exactly [8,9]. 
The numerical values of the first few coefficients (for Nf = 4) are given in Tabled Obviously, 
there is problem: the coefficients increase so fast with increasing powers of a s (rrib) (corresponding 
to subleading logarithms) that the sum in the exponent cannot be determined. The result of 
conventional Sudakov resummation with fixed logarithmic accuracy (see Section 2.1 in Ref. [8]) 
is shown in Fig. ^ It does not reach perturbative stability. 
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Table 1: The coefficients c„,fc in Eq. @ that are known exactly (NNLL accuracy). Rows (top to bottom) and 
columns (left to right) correspond to increasing powers of n (power of a a (m)) and k (power of In N), respectively. 

The three diagonals are LL, NLL and NNLL, respectively. 
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1096.06 515.20 166.25 34.89 9.04 

20399.23 10078.43 3231.40 793.25 131.33 25.95 
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Table 2: The coefficients c n ^ in Eq. 10 computed in the large-/3o limit. Columns and rows are as in Tabic Q 



To understand the reason for this divergent behavior it is useful to consider the large-/?o 
limit, where all the coefficients in the exponent can be computed. The numerical values of the 
first few coefficients (for Nf = 4) in this approximation are given in Tabled We observe that: 
(a) the result is similar to the full QCD matrix of Table ^ (b) further subleading logarithms 
have yet larger coefficients. We therefore deduce that the divergent behavior in Table^is related 
with running-coupling effects. This is a general phenomenon [10-14]: the large-order behavior 
of Sudakov exponents is dominated by infrared renormalons. The increasing coefficients of 
subleading logarithms build up factorial growth that renders the series non-summable. Since in 
the case under consideration the divergent behavior hits already at low orders it seems that any 
attempt to use this expansion is doomed to fail. 
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Figure 1: The moment-space Sudakov factor Sud(iV, mt), and the corresponding differential spectrum, where 
the exponent is computed to LL, NLL and NNLL accuracy , as dotdashes, dashes and dots, respectively. In each 
case the spectrum is matched to the full NLO result for the magnetic interaction O7. The curves end at N — 33 
and at x ~ 0.97, respectively, where the Landau singularity in gi(X) appears. 



3 Taming the divergence by Dressed Gluon Exponentiation (DGE) 

It was recently shown [8] that, despite the divergent behavior described above, resummed per- 
turbation theory a la DGE does lead to quantitative predictions for the B — ► X s j spectrum: 

dTjE^) = Mb f c+ico dN /2E 7 \ ~ N m£_ dN PTjPV / 2EA ~ N 
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where the first inverse— Mellin integral is written in terms of physical spectral moments: 

Fn = I dE ~< (jk) N1 J- = ^ e " (w)An/, ' j ^ N - !) a / m ^)' ( 5 ) 

with A PV = Mb — Tn^, and r^ r ' PV corresponds to the perturbative moments of Eq. ©, where 
the renormalon ambiguity is regularized using Principal Value (PV) Borel summation. In the 
second inverse-Mellin integral in Eq. (jlj) we assumed T ~ 1, an approximation that is valid 
away from the endpoint. Computing the spectrum by resummed perturbation theory involves: 

(a) Defining the Sudakov exponent in T^f (as an analytic function of N in the complex N 
plane) using PV Borel summation: 
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where T{u) = (u5) uS e~ uS /T(l + u5) with 5 = f3x/Po an d where the Borel transforms of the 
quark distribution (soft) and the jet anomalous dimensions, Bs(u) and Bj(u), respectively, 
are defined in Section 2.2 in Ref. [8]. 

(b) Applying the same prescription to the quark pole mass m& in Eq. (jlj) — or, equivalently 
to A in Eq. (J5J) — when computing it from the short distance mass (e.g. from m^ s ). 

(c) Performing an inverse-Mellin integral according to Eq. ©. 



In general, the calculation of the Sudakov exponent by DGE differ from conventional Su- 
dakov resummation (truncation of the sum over / in Eq. ©) in two respects. First, on the 
purely perturbative level, by using Eq. © and incorporating the exact analytic results for the 
anomalous dimensions in the large-/3o limit [14], 

B s (u) = el u (l-u) + 0(1//%), 
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the resummation of running-coupling contributions is performed to all orders (rather than to 
a given logarithmic accuracy) making the result renormalization-group invariant [10]. The 
power-like infrared sensitivity, which shows up in Eq. © through the divergence of the series, 
explicitly appears in Eq. (jHJ) as Borel singularities at integer and half integer values of u. This 
brings us to the second difference, which is an important advantage of DGE, namely having a 
definite prescription to make power-like separation between perturbative and non-perturbative 
contributions. The PV prescription we choose guarantees that the perturbative moments are 
real valued: Sud(iV, m&)| PV = [Sud(iV*, m&)| PV ]*. This choice also defines the non-perturbative 
power corrections. 

Another key property of this calculation is the exact cancellation [14] of the leading, it = 1/2 
infrared renormalon ambiguity b between the Sudakov factor and the quark pole mass in Eq. (|4"|). 
This means that the infrared sensitivity, which severely limits the accuracy of the calculation 
in a conventional perturbative approach (Fig. HJ, is nothing else but an artifact of perturbation 
theory — specifically, of using an on-shell heavy quark in the initial state — and it is fully 
removed if renormalon resummation is applied in both the Sudakov factor and the pole mass. 

To understand this cancellation consider the process-independent definition of the quark 
distribution in the meson as the Fourier transform of (^B(Pb) [^(y)7 + $> y (0,y) ^(O)]^ B(Pb)^ 

with respect to y~ where y is a lightlike vector, & y (0,y) is a Wilson line in this direction, 
and [i is a (dimensional-regularization) ultraviolet renormalization scale. Let us compare this 
distribution to its perturbative counterpart, namely the quark distribution in an on-shell heavy 
quark defined by the Fourier transform of (b(pb) [^(y)7 + $y(0, y) ^(0)]^ &(p&)V While the 
former is unambiguous, the latter involves an on-shell heavy quark external state, which is 
unphysical. While the moments of the latter distribution are infrared- and collinear-safe as 
far as logarithms are concerned, their infrared sensitivity renders them ambiguous at power 
accuracy. The large— iV limit corresponds to the replacements iP^y~ — > N and ip^y~ — ► N 
in the two matrix elements, respectively. One then obtains the following relation between the 
two soft functions [14]: 

S(M B /N; y) = S PT (m b /N; y) e -^)K/M B jr {{N _ 1)A /M B ), (8) 

" v ' 

leading renormalon cancels out 

where power corrections on the soft scale are split into two categories: kinematic ones, which 
involve A = Ms—uib an d therefore strongly depend on the quark-mass definition, and dynamical 
ones (J 7 ), which describe the Fermi motion of the quark in the meson. The OPE shows [1,2] 
that the latter begin at o(((N — 1)A/Ais) 2 ) . As indicated in Eq. (jHJ) the cancellation of the 
leading renormalon ambiguity involves only the kinematic term and therefore it can be realized 
without using any non-perturbative input on the quark distribution in the meson. An explicit 
calculation of the quark distribution in an on-shell heavy quark in the large-/?o limit c shows that 

6 This cancellation concerns all the moments and it comes on top of the cancellation of the u = 1/2 renormalon 
in the total rate [18]. 

c Such calculations were performed in a process-specific context [14] for B — ► X s ^ and B — > X U W as well as 
using the process-independent definition of the quark distribution function [9] . 




E 7 (GeV) 



Figure 2: The moment -space Sudakov factor Sud(iV, mj) (left), and the corresponding differential spectrum in 
b — > X s *y (right) as computed by DGE with perturbative expansions of NLL (dotdashed blue line) and NNLL 
(full red line) accuracy. In the latter also the residue at u — 1/2 is fixed. The DGE spectra are both matched 
with the full NLO. For comparison the conventional NNLL result is also shown (dashed green line). 

the leading renormalon ambiguity in the exponent in S PT (m,b/N; /i) is indeed equal in magnitude 
and opposite in sign to the one in the pole mass [15, 16], confirming the cancelation in Eq. jSJ). 

Returning to the observable, the replacement of the perturbative soft function S PT (mb/N; fj.) 
in Eq. (J2J) by S(Mb/N; fj,) of Eq. (JHJ), corresponding to the quark distribution is the meson, yields 
the physical spectral moments defined in Eq. ©: 

T^v = H(m b ) Sud(N,m b ) e -(N-i)A/M B ^ N _ i)A/M B ). (9) 

V v ' 

leading renormalon cancels out 

These moments are free of the leading renormalon ambiguity. Given our limited knowledge of 
dynamic power corrections, as a first approximation it is sensible to neglect them altogether by 
setting T ~ 1. This leads to the second inverse-Mellin formula in Eq. @, which is nothing but 
resummed perturbation theory. 

In practice, the application of DGE requires evaluating the integral in Eq. @ despite in- 
complete knowledge of the anomalous dimensions of both the soft and the jet functions. What 
we currently know about the exponent includes [8]: 

• The analytic expressions in the large~/?o limit [14] for Bg(u) and Bj{u) given by Eq. (J7J). 

• NNLO (i.e. 0(u 2 )) perturbative expressions in QCD [8, 14] for both B s {u) [9, 17] and 
Bj{u) [13]. 

• The normalization of the leading renormalon corresponding to the value of B$(u) at 
u = 1/2. This value was determined [8], using the cancellation, based on the u = 1/2 
renormalon residue in the pole mass, which has been computed [8, 19] with a few percent 
accuracy from the NNLO perturbative expansion of the ratio between the pole mass and 
m^ s , taking into account the known structure of the singularity [20] . 

This information is sufficient [8] to constrain the moments up to N ~ 20 fairly well. This is 
demonstrated in Fig. [21 The most important source of remaining uncertainty arises from the 
unknown behavior of Bg(u) away from the origin, and in particular, for u > 1/2. The stability 
reflected in Fig. [2 stands in sharp contrast with the situation we encountered at fixed logarithmic 
accuracy (Fig. where the leading renormalon is untamed. 

Looking at Fig.Elone immediately observes that the support properties have changed. While 
at any order in perturbation theory the result is identically zero for > upon taking the 
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Figure 3: The partial B — > X 3 *y BF with _B 7 > Eq as computed by DGE compared with new data from BaBar. 
Inner and total error bars on the data points show systematic and statistical plus systematic errors (added linearly), 
respectively. The three theory curves shown (indicated on the plot by the corresponding values for A) correspond 
to different input values for m^ s that roughly cover the uncertainty range on this parameter. Note that the overall 
normalization of the theoretical BF, which is fixed here to its central value of BF(iJ 7 > Eo = ms/20) = 3.73- 1CF 4 , 

has ~ 10% uncertainty. 



PV in Eq. JSJ) the distribution becomes non-zero there. The J5~ = nib/2 boundary is a direct 
consequence of the assumption that the initial b quark is on-shell. This inherent limitation of 
perturbation theory is removed upon defining the sum of the series in moment space by PV. 
Remarkably, the PV result tends to zero near the physical endpoint. This, together with the 
observed stability, suggests that in this formulation the dynamical power corrections associated 
with the Fermi motion (J~), which were neglected here, are indeed small. 

4 Comparison with data 

Experimentally-favored observables are the partial branching fraction (BF) for E~ > Eo, 

bf ^-T^ 7 ^ ; r( El >E )=[ dE 7 ^pl, (io) 

1 B >anything J Eq arj^ 

where the cut value E is above 1.8 GeV, as well as the first few spectral moments defined over 
similarly limited range of photon energies: 



1 f JE1 dF(E 1} 



E ih>-o ~ E V h>*o = T{Ej > e )J Eo dE ^ ~1e^ i^M " E V ■ (12) 

Measurements of the average energy defined in Eq. (|11|). and the variance, re = 2 in Eq. (|12() . 
with E = 1.815 GeV were published [21] by Belle in 2004. 

The method presented in the previous section facilitates theoretical calculation of the B — > 
X s j spectrum, and thus also evaluation of the cut moments defined above, by purely perturbative 
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Figure 4: The first two cut moments of the B — > X s ^/ spectrum: (E^) (left) and {{(E^) — -E7) 2 ) (right), as a 

function of the minimum photon energy cut Eo, as calculated by DGE (varying m^ s and a a within their error 
ranges), compared with data from Belle and BaBar. Inner and total error bars show systematic and statistical 

plus systematic errors (added linearly), respectively. 

means, so long as power corrections (J-) can be neglected. Since these power corrections are 
controlled by the soft scale rrib/N, they become increasingly important with increasing Eq. 
In order to test the theory and, eventually, determine these power corrections it is therefore 
useful [8] to compare the theoretical calculation with data as a function of Eq. Such comparison 
was done for the first time during the Moriond meeting, following the talk by J.J. Walsh [22] 
who presented new (and preliminary) results from Babar; see Figs. 01 and 0J 

5 Conclusions 

We presented a new approach to inclusive decay spectra in QCD and its application to the 
B — > X s ^f spectrum. Our method, Dressed Gluon Exponentiation, incorporates renormalon 
resummation as well as Sudakov resummation and uses the PV prescription to separate between 
computable perturbative contributions and non-perturbative power corrections, which become 
important near the endpoint. 

We have shown that, despite the apparent infrared sensitivity which shows up as divergence 
of the perturbative expansion in Eq. (j3J), the predictive power of resummed perturbation theory 
is high. This is reflected, for example, in the stability of the DGE result shown in Fig. [21 
A key property underlying this predictive power is the exact cancellation of the leading infrared 
renormalon ambiguity [14]. This cancellation can be understood as a manifestation of the 
infrared- and collinear-safety of the physical spectral moments at power accuracy. 

In the calculation of the spectrum by Eq. © the leading renormalon cancellation is realized 
by using the same prescription when defining the sum in the Sudakov exponent and when 
computing the quark pole mass which sets the overall energy scale. Renormalon resummation in 
both these quantities is of course imperative for this cancellation to take place. The final result 
for the spectrum is prescription independent as far as this leading renormalon is concerned. 
Having neglected higher power corrections on the soft scale, our result does depend on using 
the PV prescription for higher renormalon singularities. While we do not assign to the chosen 
regularization any physical meaning, we find the PV prescription particularly useful. Specifically, 
the PV-regularized spectrum shares the following key properties with the physical spectrum: 
(a) the moments are real valued, (b) the resummed spectrum smoothly extends beyond the 
perturbative endpoint and, quite remarkably, tends to zero for E^ = (mb + O(A)) /2, close to 
Mb/2. This suggest that in this formulation power corrections are indeed small. 



Finally, as shown in Figs.|3]andEJ there is good agreement between the theoretical predictions 
and data from the B factories for both the partial BF and the first two spectral moments over 
the entire range of cuts where data is available. Fits can readily be performed to measure 
m™ s and to put bounds on the dynamical power corrections, which were so far neglected in 
this approach. Finally, this successful comparison with data shows that prospects are good for 
precision determination of V u b from charmless semileptonic decays using DGE. 
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